Deterministic growth model of Laplacian charged particle aggregates 
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The results of the computer simulation of the aggregates growth of the similarly charged particles 
in the framework of deterministic Laplacian growth model on a square lattice are presented. Cluster 
growth is controlled by three parameters p, E, A, where p - Laplacian growth parameter, E - energy 
of a particle sticking to a cluster, A - the screening length of electrostatic interactions. The phase di- 
agram of cluster growth is built in the co-ordinates E, A. The zones of different cluster morphology 
are selected: I-the zone of finite X-like structures, II-the zone of infinite ramified structures, con- 
trolled by electrostatic interactions, III-the zone of infinite structures with electrostatic interactions 
effectively switched off. Simple electrostatic estimations of the locations of the zone boundaries are 
presented. It is shown that in general case within the zone II the continuous change of D/, controlled 
by parameters p, E, A, takes place. In the degeneration limit when the given model transforms into 
deterministic version of the Eden model (at p = 0), the crossover from linear (Df — 1) to compact 
(Df —2) structures is observed when passing through the boundary between the zones I and II. 

PACS number (s): 02.70.Lq, 68.70.+W, 75.40M 



A number of models have been developed to describe 
the growth of aggregates with complex, ramified struc- 
ture. The most important among them are the mod- 
els of diffusion-limited aggregation fractal growth (DLA 
model), compact growth (Eden model), ballistic growth, 
dendrite, and dense-branching structure growth [Q. On 
changing the local conditions controlling the process of 
particle's sticking to a cluster, it is possible to observe the 
structures with different morphologies ||-f|. The study 
of aggregation, which takes into account the specific in- 
teractions between the particles, e.g., the interactions of 
electric or magnetic origin, is of a special interest. Such 
interactions may effectively manifest themselves in col- 
loids in electrochemical deposition phenomena [[|, 
in discharge phenomena (in particular, in ball lightning 
development from unipolar plasma ||). 

This work presents the results of the growth study of 
aggregates formed by the likely charged particles within 
the framework of deterministic Laplacian growth model. 
Such model is widely used to describe the growth of 
dendrite-like structures PJlC|]. The model of such type 
allows to obtain the deterministic analogue of DLA-type 
fractal patterns. The choice of deterministic model al- 
lows also to enhance the lattice-imposed anisotropy of 
the patterns and to suppress the noise effects. 

Our model takes into account the screening effects by 
limiting the radius of electrostatic repulsion between the 
particles as r < A. Each simulation step included the nu- 
merical solving of Laplace equation on a lattice and the 
calculation of sticking probability / (oc to the normalized 
gradient of the potential) all over the cluster perimeter. 
Before the sticking of a new particle to the point (j) of a 
cluster, the calculation of electrostatic repulsion energy 
U(j) between the new particle and all the other parti- 



cles of cluster formed on the previous evolution step was 
carried out: 
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where q is a particle charge, r^- is a distance between the 
i particle of the aggregate and the j-point, = 1 at 
A and = otherwise. 

New particle joined the cluster at the (j)-point of 
perimeter only when two conditions were satisfied: a) 
/ > p in this point, and b) U(j) did not exceed some 
critical value E. Here, p is the Laplacian growth pa- 
rameter, < p < 1. Note that E corresponds to the 
short-range energy of single particle attraction to a clus- 
ter, and that this attraction may result from the Van der 
Waals interactions. 

The growth of clusters took place on the two- 
dimensional square lattice. The distances were measured 
in the lattice units, and all the particles had a charge 
q equal to 1. Thus, we may consider that A, E and U 
parameters introduced above are dimensionless. The lat- 
tice size comprised 1000x1000 and the number of cluster 
particles did not exceed 5 • 10 5 . 

So, in this model the three parameters, A, E,p con- 
trol the growth of clusters. In the regular determinis- 
tic Laplacian growth model, assuming no electrostatics, 
the growth of clusters is controlled only by the value of 
parameter p and their size is not limited. At p = 0, 
our model corresponds to deterministic Eden model and 
gives the compact clusters with fractal dimensionality 
Df = 2. At < p < 1, the structures with fractal di- 
mensionality develop; their dimensionality decreases with 
the increase of p ||, Jl0[ ] . For p = 0.3, which case will be 
analyzed in more details below as an example, we have 
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obtained, Df Q = 1.49. 

The switch-on of electrostatic interactions results in 
considerable complication of the growth patterns. In our 
case, both finite and infinite clusters may form, depend- 
ing on the values of \E,p parameters. 

Figure [l] presents the phase diagram in A — E coor- 
dinates, where the zones of different cluster morphology 
are distinguished. The boundary between the zones II 
and III is shown by the line 1 for p = 0.3 and by the 
line V for p = 0. The cluster patterns from the different 
zones of diagram (marked by the symbol ■ and relevant 
Latin characters) are presented on Fig.||. All clusters 
have specific anisotropic square-like shape due to deter- 
ministic version of a growth model for clusters growing 
on a square lattice and in a given presentation of clusters 
the lattice axes are directed along the square diagonals. 

In the zone I, the growth of finite X-like structures 
is observed, and in the zones II and III the growth of 
infinite clusters takes place. The maximal cluster size 
limitation observed in the zone I may be explained as fol- 
lows: at high enough A values (or, correspondingly, small 
enough E) the electrostatic repulsion energy at the points 
of new particles sticking U continuously grows with clus- 
ter growth. At some critical value of cluster size the 
realization of condition U > E occurs at all the sites of 
cluster perimeter, which means the halt of the cluster 
growth. The development of non-ramified X-like struc- 
tures in the zone I may be explained by almost complete 
control of the cluster growth by electrostatic repulsion 
forces in the given area of A, E values. Under such condi- 
tions, the particles stick to the cluster predominantly at 
points most distant from the cluster center, where the en- 
ergy of electrostatic repulsion is the lowest. In this case, 
the growth of clusters proceeds mainly along the lattice 
axes, i.e., on the lines directed along the diagonals of the 
square. 

At A = oo, only the growth of finite clusters with the 
size controlled by parameters E and p takes place, and 
the ramification of structures is also possible. The dash- 
filled area on Fig.[j] approximately corresponds to the 
zone of transition to the areas of infinite cluster growth. 
Note that the location of this zone practically does not 
depend on the value of the parameter p. 

The zone of infinite clusters, where fractal, dense- 
branching and dense structures develop, should be di- 
vided into two areas - II and III. The necessity of such 
division follows from the detailed analysis of fractal and 
morphological properties of clusters in these areas. 

The estimation of the fractal dimensionality Df of the 
deterministic growth patterns requires some caution. We 
used the "sand-box" method for Df estimation and ana- 
lyzed the number of particles in the squares of ascending 
size ]H| . Application of the squares instead of circles 
commonly used in this method allows to consider more 
precisely the effects of cluster anisotropy imposed by the 
symmetry of square lattice. Besides, as it is seen from 
Fig.^, the pattern structures may be subdivided quite 
clearly, especially at large E values, into the dense X- 



like nucleus (with Df = 2.0) and ramified peripheral 
part. We calculated the fractal dimensionality only for 
ramified part of a pattern. 

Figure || presents dependencies of Df on E for A = 
10, 100 (curves 1 and 2) and D f on A for E = 100 (curve 
3) in the case of p = 0.3. When moving across the phase 
diagram (Figjl]) along the E axis at some constant value 
of A (same as at A = oo), the clusters with fractal di- 
mensionality Df = 1 are observed in the zone I. The 
Df sharply increases on passing into the zone II; further 
on, the decreasing of Df takes place and at some point 
E = En -in, which corresponds to the point of transi- 
tion into zone III, the fractal dimensionality reaches the 
value Df = 1.49. As it was mentioned above, this value 
corresponds to the fractal dimensionality of clusters that 
grow at p = 0.3 without electrostatic interactions. When 
moving along the A axis at some fixed value of E, the 
clusters characterized by the constant value of fractal di- 
mensionality, Df = 1.49 are found in the zone III. But at 
some point A = \n-ni, which corresponds to the point 
of transition into the zone II, a rather sharp increase of 
Df takes place. 

The similar behavior of the fractal dimensionality is 
observed for other non-zero values of p, though localiza- 
tion of the boundary between the zones II and III de- 
pends on the p value. The Df behavior somewhat differs 
at p = 0. In this case, the clusters of the zone I are char- 
acterized by Df = 1 (same as for A = oo), but Df = 2 
in the zones II and III. Thus, on passing from the zone I 
into the zone II across the dash-filled zone, a crossover be- 
tween the linear (Df = 1) and dense structures (Df = 2) 
is observed. The clusters of the zone II have a dense- 
branching morphology and their structure becomes more 
compact as far as the system approaches the boundary 
of the zone III (the curve V on the Fig.Q). Zone III is 
characterized by the growth of compact clusters of the 
square shape only, which is practically equivalent to the 
cluster growth conditions with switched-off electrostatic 
interactions. 

As it follows from the above discussion, zone II of the 
phase diagram corresponds to the strong effect of electro- 
static interactions on a cluster growth, whereas zone III is 
the zone of cluster growth under completely switched-off 
electrostatic interactions. Here, we propose the simple 
electrostatic estimations for explanation of the observed 
behavior. 

The absence of any considerable effect of electrostatic 
interactions on the growth of clusters inside the zone 
III may occur only if the points of cluster-growth halt, 
complying with the condition of U > E, are absent on 
the cluster perimeter. Then, the boundary between the 
zones II and III may be estimated out of the condition 
Eii-in < U . The highest electrostatic repulsion is pecu- 
liar to perimeter points coinciding with the middle points 
of the square sides (Figj2|), since these perimeter points 
are the nearest to the cluster center. To calculate the 
electrostatic repulsion energy U in these points, we pass 
from summation to integration in relation dl|) and assume 
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the ideal density profile for fractal clusters on receding 
from the center 



D f 



p(R) w R 



where R is the distance from the cluster center to the 
middle-point of the cluster side. 
The calculations show that 



E II - III KiU = n\ D '- 1 f(D f ), 



where f(Df) is the function of cluster fractal dimension- 
ality, D f . Figure ^ demonstrates the form of this func- 
tion f(Df), obtained through the numerical integration. 

Equation (||) defines the theoretical boundary between 
the zones II and III., This boundary is shown on the 
Fig|] for the case of p = 0.3 by the dashed line 2. We see 
that this theoretical boundary line goes above the real 
boundary shown by the line 1, which is not surprising 
since the density profile assumption @ does not take 
into account the existence of the compact nucleus. 

On p = 0, i.e., for the case of deterministic analogue 
of the Eden model, the growth of compact clusters is 
observed, and we obtain from (0) 
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Condition (Q) defines the boundary between the zones 
of compact clusters growth and dense-branching mor- 
phology. In this case, the theoretical boundary defined 
by (Q) perfectly describes the simulation data (line 1'). 

Authors thank Dr. N. S. Pivovarova for valuable dis- 
cussions of the manuscript and help with its preparation. 
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FIG. 1. Phase diagram of Laplacian charged particle ag- 
gregates for deterministic growth model presented in the fol- 
lowing co-ordinates: screening length, A against sticking en- 
ergy, E. The zones of different pattern growth morphology 
are distinguished: I is the zone of finite X-like structures, II 
is the zone of infinite ramified structures with cluster growth 
controlled by electrostatic interactions, III is the zone of in- 
finite structures with effectively switched-off electrostatic in- 
teractions. Symbol ■ marks the points corresponding to the 
clusters presented on the Fig. g, which were grown under 
the following conditions: A = 100 — E — 14(a), 40(6), 100(c); 
A = oo - E = 50(a'), 100(6'), 500(c'). The solid lines 1 and V 
show the boundaries between the zones II and III for p=0.3 
and p=0, correspondingly. The dashed line 2 shows the the- 
oretically predicted position of this boundary for p=0.3. 



FIG. 2. Typical cluster patterns for the model of deter- 
ministic Laplacian growth of charged particles aggregates ob- 
tained at different values of p,E,\ (which are shown on the 
figure). Characters a, a', b, b', c and c' correspond to the 
points, shown on the phase diagram (Fig.pl). 



FIG. 3. Dependencies of the fractal dimensionality Df 
versus the sticking energy E at the fixed values of 
A = 10(1), 100(2) and versus the screening length A at the 
fixed value of £7=10 (3). 



FIG. 4. Profile of f(D f ) from the equation (3). This func- 
tion is used for theoretical estimation of the boundary be- 
tween the zones II and III on the phase diagram (Figjl]). 
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